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Abstract 

Rigorous results on the nonlinear dynamics of a dilute Bose gas with a neg- 
ative scattering length in an harmonic magnetic trap are presented and suffi- 
cient conditions for the collapse of the system are formulated. By using the 
virial theorem for the Gross-Pitaevskii equation in an external field we ana- 
lyze the temporal evolution of the mean square radius I = f r 2 \ ip | 2 dV of 
the gas cloud. In the 2D case this quantity undergoes harmonic oscillations 
with frequency 2u;o- It implies that for a negative value of the energy of the 
system, the gas cloud will collapse after a finite time interval. For positive en- 
ergy the cloud collapses if the initial conditions correspond to a large enough 
amplitude of oscillations. Stable oscillations with a small amplitude are pos- 
sible. In the 3D case the system also collapses after a finite time for a state 
with negative energy. A stringent condition for the collapse is also derived. 
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The recent observation of Bose-Einstein condensation of atoms of 87 Rb [|TJ, 7 Li and 
23 Na confined in magnetic traps has opened a new important field of investigation of 
quantum phenomena. In particular, a new feature of these systems is that they are sig- 
nificantly inhomogeneous so one can directly observe effects of quantum uncertainty on a 
macroscopical scale. The formal evidence of such a situation is that dynamics of these non- 
uniform gases at low temperatures can be properly described by a non-linear equation for 
the condensate wave function i/j(r,t) Q, f|: 

•ftj^(r,t) = (-^- + U(v)+g\^(r,t) | 2 ^M). (1) 

Here U is the confining potential and g = 4n ^ a , where a is the s-wave scattering length. 
This equation is quite classical in its merit but contains explicitly the Planck's constant H . 
(We will put H = m = 1 below.) 

Experimental evidences of BEC in magnetically trapped alkali vapour stimulated recent 
numerical investigations of Eq. (|1|) in static, linear and nonlinear regimes. || - JT0[ 



The physical situation is very different for 87 Rb, 23 Na on one side and 7 Li on the other. 
In the first case a > and the homogeneous state of the gas is stable. For 7 Li the scattering 
length a is negative (a = —1.43 x 10~7 cm), which corresponds to an effective attraction, and 
the gas possesses negative compressibility. The uniform state in this case is unstable and 
the gas can undergo a "collapse" which is an unrestricted contraction. (The idea of such a 
collapse was first introduced by V.Zacharov for a formally similar system of Langmiur waves 
||1 1|| .) Fetter however has recently shown J7j, using a variational approach, first employed 
by Baym and Pethick ||, that a system of a finite number N of atoms in a trap can be in 
a metastable state for N below some critical number N c . The same conclusion has been 
obtained by Ruprecht et al. ||. Surprisingly the number of atoms reached in the experiment 
appears significantly larger than the theoretical estimate of N c . 

In such an uncertain situation it is important to understand better the conditions and 
dynamics of the collapse. 

In this paper we will present some rigorous results relative to the behaviour of the solution 



of Eq. ([!]). Particularly we will formulate a sufficient condition for collapse in the gas. We 
will analyze the temporal evolution of the quantity / = / r 2 | ip | 2 dV, which is proportional 
to the average square of the radius of the gas cloud. This "virial theorem" has been used 
by Vlasov et al. to analyze the problem of the collapse, described by Eq. (jlj), in the absence 
of the confined field ||12|| . (See also more recent papers |Tj| - |TjJ and references therein.) 



Note first of all that Eq. ([!]) provides conservation laws for the energy E of the system: 
E = J (- | W I 2 +| I V I 4 +U I V I 2 ) dV = E + JlJ\ij\ 2 dV (2) 

and for the number of atoms N = J | if) \ 2 dV . The last conservation law is expressed by 
the continuity equation: 

The equation for the momentum density in the presence of an external field can be written 
in the form: 

dji ^ dliik \ip \ 2 ^ (A) 
dt dxk dxi 

We will not write explicitly the well known expressions for the currents ji and U.^ and 

we note only that in 2D and 3D we have an exact identities 

J Tl ti dV = 2E (5) 



and 



Ju u dV = 3E -^J | | 2 dV (6) 

respectively. 

Repeating differentiation of I with respect to t and integrating by parts one gets: 

d -L = 2j xdid V,§ = 2jlU,JV-2j X ™^f d V. (7, 
We first analyze the two-dimensional motion of the gas in an isotropic cylindrical trap 

2 2 

with a confining potential U = l ^— 1 r 2 = x 2 + y 2 . In this case Eq. (|?p takes a surprisingly 
simple form: 



fL = -A^I + AE. (8) 

whose general solution is of the obvious form: 

I = Acos(2cj £ + 7) +E/lUq. (9) 

It is worth noting that the parameters A, a and E of the solution @ can be calculated using 
the initial distribution r/>(r, 0) with the help of Eq.([|) and the first equality in Eq.fl7|). 

Eq. (H) describes harmonic pulsations of / with constant frequency 2uq in spite of the 
nonlinear character of Eq. ([!]). (It is possible to show that small axial symmetric oscillations 
of the system have the frequency 2ujq [0.) Furthermore, one can see that for a negative 
value of the energy the positive quantity / approaches zero after a finite time smaller than 
7t/2uq. This means complete contraction of the cloud. On the contrary, for a positive value 
of the energy pulsations with a small enough amplitude A < E/uq do not result in a collapse 
p~9|] . (See also the discussion in the concluding part of the paper.) 

Thus increasing the energy improves the stability of the cloud. For example an increase 
in the stability has been noted in [2y| in the presence of a vortex line. 



Equation @ has the same form both for attractive and repulsive interactions. In the 
latter case however the collapse does not occur. The confining potential energy UqI/2 is 
then always smaller than E. This inequality ensures that the collapse cannot be reached. 

One can suggest that under the experimental conditions some motion of the superfluid 
gas exists that increases the energy of the system. This means that the cloud is not in an 
equilibrium but rather in a nonstationary stable state. 

Results for the 3D case are not so conclusive. Consider an anisotropic harmonical trap 
with a confining potential of the form U = aik * iXk . One gets from ([7]) and (ffl): 

^ = 6E - 5 J a ikXi x k dV - 2K, (10) 

where K — ^ J \ \7ip | 2 dV is the kinetic energy. Eq. ( |l0l) gives the inequality: 

fL + 5u 2 m I-6E<0, (11) 



where is the smallest eigenvalue of the otik matrix. We can prove again that the system 
collapses after a finite time for a negative E. Consider the auxiliary equation 

g + hull -QE = 0. (12) 

For negative E the quantity / reaches, according to this equation, the value zero after a 
time smaller than . The negative term in the left hand side of Eq.([H]) accelerates this 
process for the same initial conditions. The same situation takes place for a positive energy 
for a large enough amplitude. 

One can get more rigorous sufficient conditions of the collapse using the powerful method 
developed in [I3|- [IB|. In this method an important part is played by two unequalities 
imposed on the kinetic energy K [13|, fi~5ll . The first one is the uncertainty relation between 
I and K: 

IK > ^N 2 . (13) 

The second inequality can be written as 

\ J I ip I 4 dV < PN 1 / 2 ^ 2 , (3 = 0.0575. (14) 

Combining ( |I~3"D and (|14[) one immediately gets an inequality for the energy E as a function 
of K: 

E(K) >K + —u 2 n N 2 K~ 1 -P\g\ K^ 2 N 1 / 2 . (15) 
16 

This equation demonstrates a remarkable resemblance to the expression for E(d~ 2 ) which 
has been obtained by Fetter using a wave function of the form ip = Cexp(— r 2 /2d 2 ). (See 
Eq.(^) in (J7|.) For such a function the kinetic energy is proportional to d~ 2 . It easily checked 
that E(d~ 2 ) (expressed as a function of K) is larger that the right hand side of fllSD . The 
equation ( |15|) permits us to formulate the stringent collapse conditions we are looking for. 
Introducing dimensionless variables e and k defined by 

E = -N(j m e, K = -Nu m K, (16) 



we can rewrite this equation as: 

e > k + k- 1 - 1.18-ct = f(K,a), (17) 
o 

where a is a dimensionless parameter 

a = V^iV | a | wjf (18) 

The parameter a is analogous to the parameter a defined in 0. 

Behaviour of the function f(n) depends on the value of a. For o > o c = 0.454 the curve 
is a monotonous one (the curve (a) in Fig.l). Let the point e* , k* to be an intersection of 
the curve e = /(/«) and the straight line e = k/3. Let us to be 

E < E* = -Nu m e* > 0. (19) 

Then the kinetic energy K > K* = 3E* for any t. This means that equation ( p!U|) leads to 
the inequality 

*L + 5^/ - 6{E -E*)< 0, (20) 

This formula obviously demonstrates that at condition (ITS) we get the collapse. 

We shall now discuss the case where a < a c . The curve e = /(/«) has in this case a 
minimum e m and a maximum e« (the curve (b) in Fig.l). Let us consider an initial state, 
where 

e m < e < e M , n x < k < k 2 . (21) 

Since any solution of (P evolves continuously, the kinetic energy k will stay in the same 
range at every time. The uncertainty relation (O) gives then 



Thus the mean square radius stays restricted and collapse of the cloud as a whole is impos- 
sible. 



If E < E* the cloud will collapse as in the previous case. The evolution of the cloud for 
E > E* and K out the interval (|2T| ) demands additional investigation. 

The method we explored is also useful in the two-dimensional case. In this case instead 
of ([13]) and (Tjj ) we have 



IK > \n 2 ,\ f | i) | 4 dV < (3NK,(3 = 0.17. (23) 

We introduce dimensionless variables e and k defined by: 

E = -Nuj m e, K = ]^Nu m K. (24) 

Then we obtain the following inequality for the energy: 

e > k{1 - N/N c ) + k- 1 = g(K } N), (25) 

where N c =l//3\g \ The shape of the right hand side of this inequality is different for 
N > N c and for N < N c . (See curves (a) and (b) in Fig.2 respectively.) In the latter 
case function g tends to +oo as k tends both to or oo. Hence the kinetic energy for a 
given value E stays restricted by some values K\ and K2 and the total collapse is impossible. 
This means that for N < N c an arbitrary initial function ip(r, 0) automatically ensures such 
initial conditions for Eq. (§) that the collapse to be absent. In the opposite case, when 
iV > N c , it is impossible to draw by this method any general conclusions. 

I would like to thank S. Stringari for various enlightening discussions and S. Giorgini for 
useful remarks. I am grateful for hospitality to the Department of Physics of the University 
of Trento where this paper has been completed. 
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FIGURES 

FIG. 1. The function /(/c) defined in (17) for a = 0.51 > a c (curve (a)) and for a = 0.29 < a c 

(curve (b)). Admited states of the system are above the bold curve. The dashed line is e = k/3. 
FIG. 2. The function g(n) defined in (25) for TV = 0.4iV c (curve (a)) and for N = 1AN C (curve 

(b)). Admited states of the system are above the corresponding curve. 



10 



